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ABSTRACT 

We present an analysis of far-infrared dust emission from diffuse cirrus clouds. 
This study is based on serendipitous observations at 160/xm at high galactic lati- 
tude with the Multiband Imaging Photometer (MIPS) onboard the Spitzer Space 
Telescope by the Spitzer Infrared Nearby Galaxies Survey (SINGS). These obser- 
vations are complemented with IRIS data at 100 and 60/im and constitute one of 
the most sensitive and unbiased samples of far infrared observations at small scale 
of diffuse interstellar clouds. Outside regions dominated by the cosmic infrared 
background fluctuations, we observe a substantial scatter in the 160/100 colors 
from cirrus emission. We compared the 160/100 color variations to 60/100 colors 
in the same fields and find a trend of decreasing 60/100 with increasing 160/100. 
This trend can not be accounted for by current dust models by changing solely 
the interstellar radiation field. It requires a significant change of dust properties 
such as grain size distribution or emissivity or a mixing of clouds in different 
physical conditions along the line of sight. These variations are important as a 
potential confusing foreground for extragalactic studies. 
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Introduction 



The InfraRed Astronomical Satellite (IRAS) showed for the first time that e xtended in- 



frare d emission was present at high galactic latitude, far from star forming regions (ILow et al. 
1984 ). In these diffuse regions, clouds are optically thin to stellar radiation and the radi- 



ation field is relatively uniform whi c h results in very limited variations of dust equilibrium 



temp erature (IBoulanger et al.l Il996l ; lArendt et al.l Il998l ; lLagache et al.l 1 19981 ; ISchlegel et al. 
19981 ). These high latitude cirrus also show a tight correlation between their infrared 



emission (100/im to 1 mm observed by DIRBE and FIRAS) and the HI column density 
( IBoulanger et al.l 119961 ) and the dust emission is well characterized with a constant dust 
emissivity per hydrogen atom {r/N H = 10~ 25 (A/250)~ 2 cm 2 ) close to the value expected 
from models of interstellar dust grains (jDraine &: Led Il984l ) . At shorter w avelength, the 
smaller dust grains emission is characterized by a ratio of Ieo/im/ hoo^m ~ 0.2 (ILaureijs et al. 



199ll ; lAbergel et al.lll996l ; IBoulanger et al.ll2000l ). All in all, dust emission from local cirrus 
is then seen as rather homogeneous and simply characterized on large scales. However, little 
is known about the dust properties (e.g. optical properties for absorption and emission, 
distribution,. . . ) in these high latitude clouds, at resolutions smaller than the DIRBE beam 
(0.7°). 

Smaller scale analysis of infrar ed colors have been d one on individual r e gions and show 



clear variations of dust properties. ILaureijs et al. I Jl996h and 



Bernard et al. 



Abergel et al.l (11994 ) observed 



(119991 ) studied the far infrared 



a decrease of Imum/ hoo^m toward dense clouds, 
emission at the arcminute scale in the Polaris flare with IRAS, ISOPHOT and PRONAOS 
(200 to 600/zm) in a region where extended emission from cirrus is detected as well as a denser 
structure. The spectrum of the extended cirrus indicates a low dust temperature associated 
with a low 60/100 fim ratio. This was also observed in the Polaris flare to ward moderately 
dense regions (Ay y 1) and in a denser filament in the Taurus complex (jCambresy et al. 



200ll ; IStepnik et al.l 120031 ). It might be explained by the formation of large dust aggregates 
through the adhesion of small dust particles onto the surface of larger grains, leading to a 
change of dust emissivity properties. In the dense regions the very small grains seem to have 
disappeared almost completely. However, all these observations were restricted to individual 
regions, most of which are much denser than the diffuse local interstellar medium seen at 
high galactic latitudes. 

By comparing near-infrared extinction and extinction deduced from far-infrared dust 
emission in the whole anticenter hemisphere, ICambresy et al.l ( 120051 ) observed a discrepancy 
between the two quantities in regions above 1 mag. This effect is also interpreted by a 
change of dust emissivity due to the presence of fluffy grains and the grain-grain coagulation 
scenario was therefore extended to larger regions. 
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Kiss et al.l ( 120061 ) analyzed the far-infrared emission properties in a large sample of 
interstellar clouds observed with ISOPHOT with respect to extinction in regions of the order 
of 100 arcmin 2 . They find variations of the far infrared dust emissivities in the coldest (12K< 
T d <14K) and densest regions that are consistent with a dust grain growth scenario. But 
they also observe changes of the dust emissivities in the warmer regions (14K< Td <17.5K) 
and interpret them as an effect of mixing along the line of sight of components with different 
temperatures or a change of the dust grain size distribution. However, a fraction of their 
sample was chosen on the basis of high brightnesses in the IRAS bands and could therefore 
be biased toward regions with enhanced small grain emission. 

Extinction measures toward high galactic latitude sightlines show a substantial fraction 
of low Ry = Ay/E(B — V) values, also indicative of enhanced relative abundances of small 
grains. However, with the lack of longer wavelength measurements at small scales, it is 
difficult to relate these variations to possible changes in the dust grain properties. 

The photometric data from the Spitzer Space Telescope enable us to have access to 
sensitive observations up to 160/im. Among the large programs, the Spitzer Infrared Nearby 
Galaxies Survey (SINGS) observed a sample of 75 nearby galaxies in photometry with the 
IRAC and MIPS instruments. The fields observed were chosen to be at high galactic lat- 
itude in order to limit the foreground cirrus contamination in the study of the targeted 
galaxy. Because the region observed was larger than the targeted galaxy, these serendip- 
itous observations are then ideal to study far infrared dust emission in a large sample 
of high galactic latitude regions. We combine these new observations with IRIS data 



( iMiville-Deschenes fc Lagachdl2005l ) at 60 and 100/im, a reprocessing of IRAS data including 
a better calibration of the infrared brightnesses and better zodiacal light subtraction. The 
goal of this paper is to study the infrared colors of diffuse local dust emission on the scale of 
a few arcminutes. 



2. The data 



The Spitzer Infrare d Nearby Galaxie s Survey fSI NGS iKennicutt et al.l 120031 ) observed 
m imagery with IRAC jFazio et al.lbooj ) and MIPS faieke et al.lhooi ) onboard Spitzer a 
sample of 75 nearby galaxies. While the IRAC images only observed the galaxy itself, a 
significant part of the MIPS observations (strips) encompass the surrounding sky. Since 
the SINGS observations were chosen to be at high galactic latitude to limit the galactic 
foreground contamination, the MIPS observations at 160/im provide a good opportunity to 
study the low surface brightness diffuse infrared emission from high galactic cirrus in a large 
number of fields at a resolution of ~ 37" . These data are complemented with IRIS data 
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( IMiville-Deschenes &: Lagachdl2005l ) at 100 and 60/un. The position of the fields on the sky 



are shown in Fig. [T] while their characteristics are summarized in Tab. [TJ 

Although SINGS observations were also done at 70/im with MIPS, the regions observed 
are offset with respect to the galaxy targeted and only a small fraction of the 70 and 160/im 
observations overlap outside the galaxy itself, making them inappropriate for our galactic cir- 
rus emission study. 24/xm observations were also available, but they are dominated by point 
sources emission as well as stronger zodiacal light. Once the point sources removed and the 
regions at low ecliptic latitude are discarded, the 24/zm brightness have a low dynamic range 
in each field and no meaningful correlation can be done with longer wavelength observations. 
This study was therefore restricted to the comparison of 60, 100 and 160/im brightnesses. 
The IRIS 25/im observations were however used together with longer wavelength in order to 
remove point sources (like galaxies) more efficiently in the observations. 



2.1. Data treatment 

The observations at 160/mi were reduced using the GeRT softwar^j] on the raw MIPS 
observations. Standard parameters were used for the reduction, but data where flashes of 
the internal source led to a significant number of saturated pixels which were removed. The 
removed data are most often positioned on the bright center of the galaxy. Other saturated 
pixels removed from processing were due to cosmic ray hits. These saturated flashes when not 
removed can bias the sensitivity of the diffuse extended emission. This step in the reduction 
may be not appropriate for the photometry of the galaxy, but significantly reduces latents 
(stripes) in the outer regions we are interested in. Each region targeted was observed twice. 
Discrepant fluxes at the same position between the two observations are removed and the 
data are combined into a mosaic for each region. 

The MIPS 160/im maps and IRIS 60/im are convolved to the IRIS 100/mi resolution 
assuming gaussian beams with FWHM of 4.0', 4.3' and 37" for IRIS 60, 100/im and MIPS 
160/mi observations respectively. 

For each MIPS strip, the galaxy and other p oint sources are detecte d in th e 25, 60 



and 100/im maps using the method described in iMiville-Deschenes et al.l (120021 ) . These 
point sources at the IRIS resolution (but with the MIPS sampling) are then smoothed by a 
gaussian kernel with a full width half maximum of 3 x 3 pixels (at the MIPS original pixel 
size) to encompass possible extended emission from these galaxies. The smoothed point 



l http : //ssc . spitzer . caltech. edu/mips/gert /index. html 
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sources are then masked in all the maps. All maps are then projected on the IRIS grid 
to avoid oversampling. The observation targeting the galaxy Holmberg IX was removed 
from the sample since the emission in the whole strip is dominated by the galaxy and its 
interaction features with nearby galaxies. We chose to remove the observation containing 
the galaxy NGC3034 (M82) which was hampered by saturation effects in the whole central 
region of the galaxy, affecting the observation globally. The observations containing the 
galaxies NGC1266, NGC2915 and M81 Dwarf B were also removed because the width of the 
region observed were too narrow to be convolved meaningfully to the IRIS resolution. We 
ended up with 70 fields of view^ at a resolution of 4.3' observed at 60, 100 and 160/mi. Due 
to uncertainties in the zodiacal light subtraction at 60/im that can dominate the flux at the 
low surface brightnesses we sample, we limited the sample at 60/xm to the 9 observations at 
high ecliptic latitude (\(3\ > 15°). 



A constant brightness of 0.78 M Jy/sr is removed fro m the IRIS 100/im maps to account 
for the cosmic infrared background (ILagache et al.ll2000i ). i.e. the emission from the distant 
unresolved galaxies (called hereafter CIB). The exact level of CIB emission has not yet been 
established at 60/im and the MIPS observations can have offsets in the calibration of the 
brightness that are not physical. To overcome the uncertainties (physical or instrumental) 
on the zero levels in the different maps, we hereafter perform the analysis of the data through 
the use of correlations (c.f. §3.ip . 

The errors on t he surface brightness are taken to b e 0.03 and 0.06 MJy/sr at 60 and 
100/im respectively (jMiville-Deschenes fc Lag ach]l2005h . At 160/zm, we take a quadratic 
combination of a constant sensitivity limit of 0.12 MJy/sr [f| and a 2% uncertainty on the 
brightness ( due to the uncertainty on the calibration factor from instrumental units to surface 
brightness (jStansberry et al.ll2007l )). 



2.2. Sample selection 

In low surface brightness regions, the variations of the infrared emission in the obser- 
vations can come from cirrus emission, fluctuations in the cosmic infrared background or 



2 Although there are 75 galaxies in the SINGS sample, some galaxies are in the same field of view: 
NGC3031 is in the same observation as M81 dwarf B and NGC5195 was observed simultaneously with 
NGC5194 

3 the sensitivity of the observations is computed for a 16s integration time per pixel using the SENS-PET 
tool, http: / /ssc. spitzer.caltech.edu/tools/senspet7l and is divided by >/~N where N — 49 is the number of 
MIPS 160//m PSF inside an IRIS PSF at 100/xm 
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from noise. Since we want to study the variations of cirrus emission only, we want to select 
observations in the SINGS sample that are dominated by the dust emission variations. 

The different contributions (cosmic infrared background, cirrus emission) to the infrared 
emission have different power spectra that can help to disentangle them. In particular, the 
cirrus power spectrum normalization depends on the mean surface brightness while the 
contribution from background galaxies does not. This dependence can be translated into 
a relationship between the mean brightness an d the standard deviation squa re, a 2 ir 



m 



a region and depends on the size of the region (IMiville-Deschenes et al.l 120071 ). For each of 
the observed field of view, we computed the standard deviation at 100/im and the mean 
brightness at 100/im (minus the average CIB contribution at this wavelength) and then plot 
the cr 2 -< B > relations hip observed at 100/xm i n our s ample. To model a cirrus , we use the 
relationship derived by IMiville-Deschenes et al.l (120071 ) below 10 MJy/sr, for a maximum 
scale length of 50' (dotted line in Fig [2]). We observe that our observations are consistent 
with the model, with a large scatter as in the original relationship. This dispersion is likely 
enhanced due to the fact that our fields of view are elongated and the size of the region 
mapped varies between field. The contribution from the CIB fluctuations can be described 
by two terms: a Poisson noise that represents the galaxies distributed homogeneously with 
respect to the resolution and a component with correlated spatial variations corresponding 
to the clustering of galaxies on lar ge scales. The c ontrib ution from the clustering of infrared 
galaxies is predicted by using the lLagache et al.l (120031 ) model for galaxy evolution, with a 
bias parameter from Lagache et al. ( 20071 ). The contri bution from the Poisson noise to the 
a 2 observed at 100/im is taken to be that measured by IMiville-Deschenes et al.l (120021 ) since 
we used the same point source detection method. However, compared to their study, we 
removed point sources applying the detection scheme at all wavelength (25, 60 and 100/im). 
This enables us to mask faint galaxies at 100/im more efficiently and the Poisson noise in our 
measurements could be lower than their measurement. Because we want to select the fields 
with the least contribution from other sources (CIB) than cirrus to the observed variations, 
this choice is therefore conservative. 

Combining all contributions (represented in Fig. [2]) to the observed variations, we deter- 
mine that a cut at 2.5 MJy/sr corresponds to o cirrus j o C ib = 1-5 (so that the total infrared 
fluctuations a tot 



(7,. 



i + a ciB are l ess thatn 20% larger than from cirrus fluctuations 
alone). The regions with a mean 100/im brightness above this threshold will therefore be 
dominated by variations of cirrus emission. In each field, we computed the mean brightness 
at 100/im as well as the standard deviation at 60, 100 and 160/xm (c.f. TabJT]). By keeping 
only the fields above the 2.5 MJy/sr cut, the sub-sample we will study in this paper is com- 
posed of 15 fields with 100 and 160/xm brightnesses, among which 9 can be studied as well 
at 60/im > 15°, see sec. El]). 
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Stellar reddenings obtained from the analysis of the Sloan Digital Sky Survey data 
enable us to put an upper limit of 1.2 mag on the extinction in these fieldfl This confirms 
that the variations in t he infrared cirrus emiss i on stu died in each region comes from diffuse 
clouds according to the Ivan Dishoeck fc Black! (119881 ) classification. 



Results 



3.1. Cirrus emission at 160 and 100/xm 

In each field of the selected sample, we plot the point to point correlation between the 
brightnesses observed at 100 and 160/im (represented in Fig [3] and H]) and apply a linear fit 
taking into account the errors at both wavelengths. This enables us to obtain for each field 
a slope corresponding to the ratio Si 60 /-Bioo unbiased by variations of the zero point level 
(residuals from the zodiacal light subtraction, absolute value of the CIB). The correlation 
coefficient and the slope derived in each region are summarized in Tab. [2j 

Large scale observations of high galactic latitude emission of cirrus with COBE were well 
characterized by a modifie d black body with a d ust temperature of 17.5K and an emissivity 
index proportionnal to v 2 (IBoulanger et al.lll996l ). Using this law, we estimate the large scale 
160/100 color for cirrus to be of -B160/-B100 = 2.0 (taking into account color corrections). This 
ratio is represented in the correlation plots (Fig. [3] and @J to guide the eye. While some 
correlations between B wo and £> 160 are in agreement with the -Bi 60 /-E>ioo — 2.0 obtained on 
large scales, clear deviations are also seen (4 fields out of 15 have a slope discrepant at a 5-6<r 
level with respect to the value of 2.0. The most extreme case is the field of NGC2976, with 
a fitted slope on the Biq versus Bi 00 correlation that is 10cr away from the 2.0 standard 
value) . 

In Fig. we compare the obtained ratios Bi 60 /B 100 to the mean surface brightness at 
100/iHi in each field (black points). We observe a large dispersion in the 160/100 colors that 
can not be explained by the error on the data or the fitting process. At 100 and 160/im, the 
interstellar emission is do minated by the emi ssion from big dust grains at thermal equilibrium 
with the radiation field (IDesert et al.lll990l ) and the Bi G0 /B 100 ratio is therefore related to 
that characteristi c dust temperature. T aking a standard emissivity of dust per hydrogen 
atom in Hi from IBoulanger et al.l (119961 ) and an emissivity index of 2, the 160/100 color 
variations we are probing can therefore be related for illustrative purposes to temperatures 



Using N{H)/A V dBohlin et alJ Il978h and B 100 /N(Hj) » 6.67 x 1CT 21 MJy/s 



this upper limit 



implies -B100 < 15- 2MJy/sr, which is fully consistent with the brightness observed in our sample 
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between 15.7 and 18. 9K for column densities ranging from Nh = 3 x 10 20 to 2 x 10 21 cm~ 2 . 
We note that these variations are consistent on average with the large scale estimate (blue 
solid line), confirming that the fields used in this study are sampling the cirrus emission 
observed on large scale. 

For a given grain size and composition, this characteristic temperature depends on the 
local radiation field strength and spectrum which depends on the presence and distance of 
nearby heating sources and on the extinction. In the framework of this model, the pres- 
ence of large variations in the 160/100/im ratio observed in our sample would suggests the 
presence of large variations in the heating of grains at small scales (variations by a factor of 
3 of the intensity of the incident radiation field). This can be surprising since at low FIR 
surface brightness and at high latitude the interstellar radiation field might be expected to 
be homogeneou s. We looked at the fa r-inf rared color temperatu re maps derived from DIRBE 
observations by lLagache et al.l (119981 ) and lSchlegel et al.l (119981 ). The regions we are studying 
appear to be reasonably representative of the high latitude cirrus given the small number 
statistics. For the sightlines covered by our sample, the FIR color variations seen in the 
DIRBE data are compatible with the variations that we observe. Our study is indeed more 
sensitive than previous observations and therefore able to probe color variations smaller than 
the uncertainties in the previous studies. 

The shape of the optical spectrum heating the grains could also affect the far-infrared 
colors: the rad iation field could become gradual ly harder with position off the g a lactic plane 
(jMattilal Il980l ). Using the cirrus model from lEfstathiou &: Rowan- Robinson! (120031 ) with 
different stellar populations heating the clouds, we checked that changes in the shape of the 
optical radiation field is unlikely to affect the 160/100 and 60/100 colors of cirrus by more 
than 20%. 

The dust equilibrium temperature depends however also on the structure of the grains. 
Grain aggregates for example cool more efficiently The temperature variations observed 
in the diffuse medium could therefore be either due to variations of the intensity of the 
interstellar radiation field or to changes in the grain structure. 

We compared our findings with different studies of far infrared emis s ion fr om the lit- 
erature: the quiescent high galactic lati tude cloud s from Idel Burgo et al.l (120031 ). the large 
sample from archival ISOPHOT data by 



Kiss et al 



(20060, the two regions in a high latitude 



cirrus MCLD 123.5+24.9 observed by [Bernard et al.l (119991) . and the quiescent filament in 
the Taurus molecular complex from Stepnik et al. ( 20031 ). Because other observations were 
obtained with different instruments, we have to interpolate the brightnesses at various far 
infrared wavelengths to estimates at 100 and 160/im. To do so, we took the dust temper- 
atures determined in each study with a brightness at 100 or 200/mi and used a modified 
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black body law with a spectral index. The power index is either taken from the study itself 
(if it was computed) or is fixed to a standard value of 2. For each region, we also compute 
the mean 100/im brightness as observed by IRIS and subtract a mean CIB contribution as 
for our observations. Despite large scatter, we observe a trend between < B WQ > and the 
160/100 color that is consistent with the idea that denser regions are colder. However, the 
effect of selection biases of these studies remains unclear. The comparison of our results with 
that from the literature (Fig. [5]) shows that previous studies in the far-infrared have been 
targeting higher Bi 60 /B W0 and < B wo >, i.e. denser and colder clouds. Due to our better 
sensitivity, our observations fill the gap at low B wo and low Biq /Bi 00 . 

For the first time, we observe interstellar dust emission at low surface brightness in an 
unbiased way (in the observing strategy) with a high sensitivity. These observations show 
that former studies on dust properties at FIR wavelength at small scale, have been biased 
toward colder and denser clouds. Our study shows that the 160/100 brightness ratios of 
high galactic cirrus clouds at small scales are consistent on average with the observations 
on large scales. However, these 160/100 colors show a wide dispersion that could arise from 
variations in the heating of the clouds or from change of the dust grain structure. In order 
to investigate further the origin of the 160/100 variations, we extend the comparison to the 
60/im data. 



3.2. Comparison to the 60/im data 



To investigate the origin of the 160/100 color variations observed in the diffuse cirrus, 
we compare the 160 and 100/zm data to the 60/im emission. The sample for this part of 
the study is however reduced to fields with an ecliptic latitude above 15° in order to avoid 
artifacts due to the uncertainties in the zodiacal light subtraction in the IRIS data. For each 
field, we determine a 60/100 color by using the same correlation technique as above. The 
correlations in each region are shown in Fig. [6] and the obtained B m /B wo are summarized 
in Tab. [3 

Fig. [7] presents the Bq /B wo ratio obtained in each field with respect to the Biq Q /B W0 
ratio. Here again, large variations are observed in the 60/100 colors that can not be ex- 
plained by the uncertainties in the data or in the analysis. As for the 160/100 colors, the 
60/100 brightness ratios are consistent on average with the "reference values" (the pink 
cross and the blue star in the figure) obtained for high latitude em ission on large scales 
(IBoulanger fc Peraultlll988l ; iBoulanger et al.lll996l ; lArendt et al.lll998l ). Furthermore, there 
is a trend of decreasing _B 60 /-Bioo with increasing Si 60 /-Bioo- 
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To try to interpret this trend , we used two mode ls of the dust grain emission at different 
interstellar radiation fields: the Draine fc Li ( 2007) model for th e Milky Way and the 



'DUSTEM" model (updated model based on iDesert et al.1 (ll990l )). The models take into 



account the shape of the IRAS and MlPS/Spitzer filters, the color corrections. For both 
models, the abundances of different grain types are kept constant. The tracks obtained are 
compared to the data in Fig. [7J The comparison shows that, if the variations of colors are 
due to variations in heating of the grains, this would imply large changes in the interstellar 
radiation field at high galactic latitude (from U ~ 0.3 to 1). Furthermore, the trend observed 
between the 60/100 and 160/100 colors is not well reproduced with the current mo dels by 



changing the interstellar radiation field alone. In that respect, the lDraine &: Lil (120071 ) model 
is however closer to the observed trend than the DUSTEM model (only 3 fields are more 
than 3cx deviant from the expected curve), but for B im /Bi 00 < 2.1 all observations show 
systematically higher B 60 /B 100 than predicted by both models, while for Bi G0 /B W0 > 2.1 all 
data points have systematically lower Bqq/B W q values than expected from both models. 

Current dust models might be missing an additional dust grain type. Such an addition 
might reproduce all color variations while changing the interstellar radiation field alone. 
Another way to interpret the observed trend is that the variations in the dust emission 
spectrum reflect spatial changes in the grain properties - composition, structure or size 
distribution. 

The equilibrium temperatures of dust grains is expected to decrease for increasing grain 
sizes and small grains (< 0.01/xm) undergo temperature excursions following single-photon 
heating that enhances the 60/zm emission. Thus regions with fewer small grains may have 
lower 60/100 ratios. The observed trend between the 60/100 and 160/100 infrared colors 
could be reproduced by changing the dust grains size distribution or composition. For exam- 
ple, enhancing the amount of small grains in regions with higher interstellar radiation field 
(i.e. higher temperatures) and reducing it at low equilibrium temperatures could reproduce 
the observed variations. 

In the same way, regions with enhanced populations of large grains may have increased 
160/100 ratios. In that case, reproducing the observations could be obtained with only 
modest variations in the starlight heating rate and shifts in the grain size distribution (fewer 
small grains and increased sizes for the larger grains at low temperatures, more small grains 
and smaller sizes for the big grains at higher dust temperatures). 

The 60/100 colors we observe therefore suggest changes of the dust properties (dust 



5 We took the model with a PAH fraction qpAH = 4.58% but we checked that this parameter does not 
influence significantly the results of this study 
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size distribution and/or composition) from one region to the other. These changes are 
related to variations in the 160/100 brightness ratio. Whether the 160/100 color variations 
require a change of the starlight intensity heating the clouds or result from the change of 
dust properties alone is unclear. The interpretation of the color variations and of the trend 
between the 160/100 and 60/100 colors is discussed further in the next section. 



4. Discussion 



Despite its rather constant color distribution on large scale, the far infrared emission 
from diffuse cirrus at high galactic latitude is observed to host large color variations on small 
scales. These variations seem related to each other (the 60/100 color decreases as the 160/100 
color increases). In this section, we will first check that these variations come indeed from 
cirrus emission and are not related to the galaxies targeted with the observations. Second, 
we will discuss possible interpretations for these large color variations and the trend between 
them. 



4.1. Extended disks in galaxies 

Because the MIPS observations were taken to observe nearby galaxies, it is legitimate 
to ask whether the infrared emission that we observe could be associated with these targets. 
In particular, Hi observations have shown that gaseous disks can extend much farther than 
the optical diameters. Dust grains could be present in these outer parts of the galaxies and 
bias our measurements. 

To avoid this extended emission from the galaxies, we were careful in masking regions 
larger than the detected emission (c.f. Sec. I2.ip . Some of the galaxies in the observations 
used in this study have been observe d in Hi observations through The HI Nearby Galaxies 



Survey (THINGS. IWalter et al.ll2005l ). We checked that the Hi diameters reported for these 
galaxies are smaller than the region masked for our study. We are therefore confident that 
the variations observed in the infrared emission between fields do not come from the targeted 
galaxies, but rather from diffuse cirrus emission. 



4.2. Physical conditions in cirrus clouds 



A possibility to interpret the variations of far-infrared colors at high galactic latitude 
is that we are sampling clouds in different physical conditions and/or composition (different 
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heating of the clouds, different dust size distribution, . . . ). The color changes would then be 
due to mixing along the sightline of these different components. 



An unbiased survey of H 2 absorption in high galactic latitude clouds by FUSE (iGillmon et al 



20061 ) has been interpreted as showing that some clouds have been compressed. The dynam- 
ical history leading to this compression may involve shock waves or strong turbulence, which 
could also lead to changes in the grain size distribution by shattering in grain-grain colli- 
sions, possibly explaining the regions of higher than average 60/100 and lower than average 
160/100 colors. 

One tantalizing possibility, in terms of mixing, is the presence of Intermediate or high 
velocity clouds (IVCs and HVCs) along the sightline. These cloud falling onto our galaxy 



could have very different dust properties (e.g. iMiville-Deschenes et al.l 120051 ) and would bias 



the measured infrared emission from more local cirrus clouds. We checked for the p resence 
of intermediate velocity clouds in the LAB Hi survey spectra (IKalberla et al.l 120051 ) in the 
direction of the fields in our sample. For about half of the sample, there is a intermediate 
velocity component seen in Hi in the sightline. Two sightlines also have a high velocity 
component. However, no conclusive trend between the fraction of the Hi in the IVCs and/or 
HVCs and the infrared colors could be seen. This may be due to the lack of resolution of the 
Hi observations (~ 0.6°), to the difficulty to disentangle IVCs and the Milky Way in some 
regions or to the small size of our sample. 



4.3. Variations in grain properties 



The grain size distribution is the result of processes such as sputtering, shattering and 
coagulation, and sightline-to-sightline variations in the wavelength-dependence of optical and 
ultraviolet extinction toward stars have already demonstrated regional variations in the grain 
properties. Whether the observed variations in emission can be fully explained by variations 
in the size-distribution alone, or whether other properties (e.g., composition or porosity) are 
also involved is uncertain. 



In denser clouds, variations of infrared colors (jStepnik et al.ll2003l ; iKiss et al.ll2006l ) have 
been interpreted with a grain coagulation scenario, combining changes of the size distribution 
with changes of the dust emissivity properties. The trend observed between the 60/100 and 
160/100 colors would be consistent with this idea. In this scenario, most of the 160/100 color 
variations would then be due to the change of emissivity of dust grains (due to changes of 
their structure), while the 60/100 colors would change with the incorporation or release of 
small grains in large dust aggregates. We could be witnessing variations of dust properties 
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due to variations in turbulent motions in the diffuse interstellar medium. Alternatively, 
dust grains in the diffuse medium could retain for some time the aggregate structure they 
had previously acquired in denser regions. So we could be seeing a sequence of regions 
corresponding to increasing time since their release from high density environments. 

Such changes in the dust size distribution and structure of grains would imply related 
variations of the UV-optical extinction curves at high galactic latitude. An extinction and 



redde ning study of stars at high galactic latitude behind translucent clouds (ILarson fc Whittet 



20051 ) shows variations in the extinction curves obtained with respect to the average curve 
for the diffuse interstellar medium. In particular, 48% of their sightlines have Ry < 2.8, 
much lower that the diffuse ISM average of 3.05. Such values are indicative of enhanced 
abundances of small grains, and these regions could have a higher that average 60/100 color. 
To test if the high 60/100 colors and low Ry values are connected, we computed the 60/100 
colors from IR IS data in a simila r fashi on to this study for a ~ 1 x 1" region around each 



sightline of the ILarson fc Whittetl (120051 ) sample. We do not observe any correlation however 
between the Ry and 60/100 color, nor between Ay and the 60/100 color. Unfortunately, no 
longer wavelength observations exist for these regions and it would be important to deter- 
mine the 160/100 colors in these regions as well and study their dependancy with extinction 
properties. This result is however a concern for the coagulation scenario as an interpretation 
of the 60/100 color variations we observe. 

Interpreting the trend observed between the 60/100 and 160/100 colors with a change 
of dust optical properties and dust size distribution remains hypothetical without further 
observations. In particular, Hi observations of these diffuse sightlines, with a high resolution 
(at least similar to the IRIS one), will be needed to determine the emissivity of dust per 
hydrogen atom and test if variations are observed and correlated with the infrared colors. 

It is presently not possible to interpret further the far-infrared color variations in terms 
of physical condition changes, grain size distribution, grain properties, etc. A larger number 
of observations of high latitude cirrus could help to probe the spatial variations. Hi studies 
of these regions at high resolution would also be crucial to probe possible variations of dust 
grain emissivities, or to check whether the velocity structure of the cirrus correlates with the 
60/100/160 colors. Finally, extinction curves on sightlines where cirrus emission properties 
have been determined would be most useful to see whether extinction properties would 
correlate with the FIR colors. 
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5. Conclusion 

We performed an unbiased study of dust emission from high galactic latitude cirrus using 
serendipitous Spitzer MIPS observations at 160/im from the SINGS survey, complemented by 
IRIS data at 60 and 100/xm. After an appropriate post-reduction of the data and a removal 
of the targeted galaxy, a sub-sample is selected so that the variations of the cirrus emission 
dominate over the CIB fluctuations in each field. 

We observe 160/100 colors in our fields that are consistent on average with large scale 
studies. However, strong variations are also observed from field to field. This paper extends 
former studies on dust properties at high galactic latitude to more diffuse, fainter and warmer 
clouds. The 60/100 color is also observed to vary significantly in the sample and there is 
a trend of decreasing 60/100 with increasing 160/100 ratios. This trend is not completely 
reproduced by current models taking only into account variations of the radiation field 
strength and requires changes in the dust properties, composition or size distribution. 

The exact origin of these variations remains unknown, but the variations of the 60/100/160 
colors may reflect variations of the grains size distribution, of grain properties in addition to 
heating variations. However, we can not completely rule out the possibility that our fields 
contain emission from matter at different heights above the Galactic plane, the juxtaposition 
of multiple components in the fields could be affecting the infrared color estimates. 

All in all, we observe unexpected variations of far-infrared colors in the supposed "ho- 
mogeneous" cirrus at high galactic latitudes. These variations are not yet understood and 
further studies will be needed to test their origin. In particular, studies on a larger area of 
sky is needed to confirm the significance of these variations and their spatial distribution on 
the sky could give new clues on their origin. 

These infrared color variations will most probably be linked with variations of the in- 
frared colors at longer wavelengths. They therefore represent an important point to study 
for the Planck and Herschel missions. Such longer wavelength observations will enable us to 
determine precisely the temperature and spectral index of the dust, and their variations, in 
high latitude regions. 

We would like to thank our referee M. Rowan-Robinson for his useful and interesting 
comments and A. Efstathiou for his help with their cirrus model. We are also grateful to L. 
Cambresy for useful discussions and work on upper limit for extinction in our fields. 
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Fig. 1. — Position of the SINGS fields (red ci rcles and blue squares) overlaid on the dust 
column density maps from lSchlegel et al.l (119981 ) centered around the north galactic pole (left 
panel) and the south galactic pole (right panel). The blue squares correspond to the fields 
selected for this study (the variation in the infrared emission is dominated by the cirrus 
component). A grid of galactic coordinates is overlayed. 
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Fig. 2. — Variations of the square of the standard deviation (related to the power spectrum of 
the signal) measured in each field with the mean brightness at 100/im. The observed values 
are compared with models for the different contributions: the infrared galaxies clustering 
(dotted-dashed line), the poisson noise (dashed line) and the cirrus variation (dotted line). 
This enables us to define a cut in 100/im brightness (the vertical black line) above which the 
cirrus variations dominate over CIB fluctuations. 
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Fig. 3.— 160-100 scatterplots for all SINGS observations with < B 100 >> 2.45 MJy/sr. 
In each plot, a canonical slope of 2.0 is represented by a dashed line (corresponding to 
a temperature of 17. 5K). A linear fit is performed on the correlation and the best fit is 
represented with a solid line. The value of the slope obtained is written in the legend. 
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Fig. 5. — Variations of the 160/100 surface brightness ratio with the mean 100/im surface 
brightness. The diamonds represent the SINGS observations. The blue line denotes a typical 
ratio of 2.0. 
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Fig. 6. — 60-100 scatterplots for all SINGS observations with < Biqq >> 2.5 MJy/sr. In 
each plot, a canonical slope of 0.2 is represented by a dashed line. A linear fit is applied and 
the best fit is represented with a solid line. The value of the slope obtained is written in the 
legend. The green point and its bold uncertainties represent the value chosen to represent 
the local background in the observation. 
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Fig. 7. — Color-color plot showing the evolution of the 160/100 and 60/100 brightness ratios 
with each other for all the the regions where < B wo >> 2.5 MJy/sr. Two reference values 
from the literature representing the colors measured on large scales are overplotted. The 
observed colors are compared to the trend obtained with two dust models: the Draine et al. 
(2007) model (dark solid line) and the DUSTEM model (blue dotted line). In both case, the 
incident interstellar radiation field intensity is the only parameter varying (the proportion 
of grain types are kept constant). 
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Table 1. Characteristics of the observations: galactic and ecliptic coordinates, average 
100/xm cirrus brightnesses in each field as well as standard deviations at 60, 100 and 160yum 
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Table 1 — Continued 



field (l,b) 


(A./3) 


Area (deg 2 ) 




< Bioo > (MJy/sr) 


CT ioo CT i60 


NGC7793 ( 5.28269 , -76.5898) 


( 344.625 , -29.1746) 


0.119 


0.046 


0.95 


0.126 0.184 



Table 2. Infrared colors in the observations, i.e. the 160/100 and 60/100 brightness ratios 
obtained from a fit of the correlations in each field 



field 160/100 color 60/100 color 



NGC0337 


2.3 ±0.2 




NGC0628 


2.6 ±0.1 




NGC0855 


2.7 ±0.5 




NGC0925 


2.4 ±0.3 


0.12 ±0.03 


NGC2976 


3.0 ±0.1 


0.15 ±0.01 


NGC3521 


1.81 ±0.07 




NGC3621 


1.86 ±0.04 


0.225 ±0.005 


NGC4569 


2.0 ±0.2 


0.31 ±0.03 


NGC4594 


2.3 ±0.1 




TOL89 


1.7 ±0.1 


0.32 ±0.02 


NGC5408 


2.22 ±0.04 


0.19 ±0.01 


IC4710 


2.0 ±0.1 


0.31 ±0.01 


NGC6822 


2.23 ±0.04 




NGC6946 


2.5 ±0.1 


0.16 ± 0.006 


NGC7331 


1.75 ±0.07 


0.238 ±0.009 



